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Hysteresis energy decreased significantly as nanocrystalline NiTi shape memory alloy was under triangular 
cyclic nanoindentation loadings at high rate. Jagged curves evidenced discrete stress relaxations. With a 
large recovery state of maximum deformation in each cycle, this behavior concluded in several nucleation 
sites of phase transformation in stressed bulk. Additionally, the higher initial propagation velocity of 
interface and thermal activation volume, and higher levels of phase transition stress in subsequent cycles 
explained the mono tonic decreasing trend of dissipated energy. In contrast, the dissipated energy showed an 
opposite increasing trend during triangular cyclic loadings at a low rate and 60 sec holding time after each 
unloading stage. Due to the isothermal loading rate and the holding time, a major part of the released latent 
heat was transferred during the cyclic loading resulting in an unchanged phase transition stress. This fact 
with the reorientation phenomenon explained the monotonic increasing trend of hysteresis energy. 



Thermomechanical NiTi shape memory alloys (SMAs) have been used extensively for intricate processes in high 
tech industries due to their excellent self-controlled and simple mechanisms\ The range of applications for 
these smart materials has covered bio-, nano- and macro- technologies, e.g., coils to treat damaged vascular 
systems, artificial heart tissues, blood clot removers, porous structures for artificial joints, nano-dampers and 
composites for aircraft elements^. In all these applications, NiTi SMA is under cyclic loadings with different setups, 
magnitudes and ratesl NiTi SMA accommodates a large amount of inelastic deformation by atomic lattice rear- 
rangements called discrete forward and reverse phase transformation^. In this transition, the parent austenite phase 
(B2) and a shear-like distorted-monoclinic martensite phase (B19') revert to each other with release or absorption of 
latent heat. A shape recovery is obtained after removal of the load (superelasticity) or upon an application of heat 
(shape memory effect). This implies that the governing parameters for a thermoelastic NiTi SMA is stress, tem- 
perature and their thermal-mechanical loading histories, e.g., rates, cycles, and loading set ups^"^\ These parameters 
have been thoroughly studied in uniaxial macro-scale torsional, bending and especially tensile conditions 
Recently, developed technologies, i.e., in-situ nanoindentation machines, have provided a great capability to reveal 
the behavior of NiTi SMAs^°"^^ and other materials^^"^° under nano-scale and atomistic level compression. 

In cyclic-compressive loadings, the temperature and initial stress of each cycle and their correlation are key 
parameters because a slight change in these parameters results in a huge change in the martensite volume in 
subsequent cycles^^'^^'^^"^l Thus, due to a significant influence on the alloy's behaviour, the thermal effect cannot 
be neglected in the case of a dynamic/cyclic loading^^. In addition, recent studies have shown a different nature of 
heat transfer phenomenon in extremely small scale compared to a macro scale^^. A different paradigm of heat 
transfer was also presented during ultrafast loadings. Nevertheless, these recent advances have been overlooked 
within literature and need to be accounted for predicting the nanoscale behavior of thermomechanical smart 
materials^'^'^^. In this paper, it is theoretically and experimentally demonstrated that the thermal effect, the 
activation volume, the orientation and the velocity of interface propagation determine the performance of 
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NiTi SMA under cyclic nanoindentation loadings. A concomitant 
hypothesis is provided for the wake of phase transformation volume 
in a nanoscale bulk according to the trend of nanoindentation curves. 
A numerical simulation is also provided with its constitutive assump- 
tions to support the results. 

Results 

Fig. 1 shows the jagged indentation curves of the cyclic loading on 
NiTi SMA with a triangular set up (inset Fig. 1). With no holding 
time, a consecutive load was instantly applied on the sample after a 
loading-unloading cycle was completed. In each cycle, the load 
started from zero to a maximum 10000 jiN, which was constant 
for all cycles. The loading and unloading rates were similar, 2 X 
10^ |iNsec~\ Each cycle had a low plastic deformation, however, 
an accumulated deformation still existed at the end. A slight increase 
in the maximum indentation depth during the 10-cyclic loading was 
observed. The maximum hysteresis area was witnessed in the first 
cycle, while the last cycle indicated the lowest dissipation energy. 
Similar results were obtained in other complementary experiments 
of 30-cycUc loadings. 

To study the response of NiTi SMA to distinct cyclic loading 
rounds, 10 indentations were performed on the same location. 
Fig. 2 shows the cyclic loading curves with the setup shown in the 
inset figure. Each cyclic routine consisted of a loading and an unload- 
ing stage at a similar rate of 2 X 10^ |iNsec"\ followed by a 60-sec 
holding time at the end of the unloading stage while the tip was in a 
holding position with no contact with the sample. The next cycle was 
performed by descending the tip on the same indented spot with the 
similar three sections of the set up. These indentations were con- 
ducted ten times. 

During the cycles, the indentation machine coordinated a new 
origin for Z (the indentation direction as Z = 0) for any fresh test 
when the tip touched the spot for a subsequent indentation. As a 
result, every indentation curve started from zero Z for a following 
cycle. As such, before stabilization, the residual depth of each cycle 
could be calculated by adding up the obtained residual depth to the 
remaining residual impression of the previous cycle. After cycle 6, 
where stabilization occurred, the residual depth did not change and 
the total residual depth remained unchanged. This fact can be 
observed in Fig. 2 where the residual depth of each individual cycle 
stabilizes after cycle 6. Similar to Fig. 1 the residual depths increased 
while the cyclic loading continued towards the 10* cycle. Contrary to 
Fig. 1, the maximum nanoindentation depth decreased from the first 
to the last cycle. 
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Figure 1 | Indentation curve of 10 consecutive cycles with a 3.61 ^im 
radius tip and 10000 ^iN maximum indentation force. Inset figure: The 
triangular set up for cyclic loading on the NiTi SMA. 
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Figure 2 | Indentation curves of 10 consecutive cycles with a 3.61 \im 
radius tip and 10000 maximum indentation force on the same 
indentation spot. Inset figure: The triangular set up of cyclic loading on 
NiTi SMA with a 60 sec holding time after each unloading stage. 

The hysteresis areas of indentation curves are shown in Fig. 3. The 
dissipated energy for the data in Fig. 1 (cyclic loading with no holding 
time) decreased from —0.23 nanojoules in the first cycle and 
approached —0.1 nanojoules in the 10^^ cycle. The increasing trend 
of the hysteresis areas for data in Fig. 2 (cyclic loading with 60-sec 
holding time) started from —0.1 nanojoules in the first cycle and 
stabilized at —0.18 nanojoules in the 6* to 10^^ cycles. The hysteresis 
areas for both cyclic set ups were similar in the 3'^-4^^ cycles and 
major stability occurred after 6 indentations. The results also suggest 
that, in contrast to Fig. 1, the hysteresis area in Fig. 2 did not decline, 
rather it did significantly increase during the ten cycles. 

Thermomechanical behaviour of NiTi SMAs. In a complete round 
of a loading-unloading stage for NiTi SMA, the hysteresis stress is 
generally defined in a simple form as^^: 



(1) 



6 and e, (Tam^ ^ma are temperature, rate (either strain or loading), 
forward phase transformation stress (austenite (A) to martensite 
(M)) and reverse phase transformation stress (M to A). During the 
loading stage on an austenitic NiTi SMA, the process of phase 
transformation generally consists of an elastic deformation of 
austenite phase, a transformation from the austenite phase to the 
martensite phase at the phase transformation stress ((Tam)^ an 
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Figure 3 | Hysteresis area of each cycle for two cyclic set ups. 
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elastic deformation of martensite phase and the plastic deformation 
of martensite phase. At the unloading stage, the resultant elastically 
deformed martensite phase reverses until the stress reaches the 
reverse phase transition {a ma)- At this stress, the reverse phase 
transformation continues until all martensitic phase transforms to 
the austenite phase. The hysteresis stress (Hin Eq. 1) represents the 
net applied stress to perform the phase transformation process from 
the beginning of the forward phase transition towards the end of the 
reverse phase transition. 

In a simple form, the temperature dependent stresses in Eq. 1 have 
linear relationships with temperature through Clausius-Claperon 
equations^ ^: 

(^am{0) = (Tam{Oo) + bAMiO - Oo), ^AM>1 (2) 

(yMA{0) = aMA{Oo) + bMAiO-Oo), bMA>l (3) 

^amj ^ma are initial temperature and Clausius-Claperon coeffi- 
cients. Eqs. (2) and (3) reveal the increasing trend of forward and 
reverse phase transformation stresses of NiTi SMA in the case of a 
temperature increase. It will be discussed in other sections that a 
higher rate results in a higher temperature depending on the 
setup^^'^^. Therefore, the rate has been separated from the temper- 
ature in Eq. 1 to highlight the possibility of a similar effect in the 
absence of one another. 

Eq. 1-3 are commonly applied for tensile tests, however, they can 
be generalized for indentation tests. The corresponding relations 
between a compression and a tensile test can be obtained when the 
normal pressure on the projected indentation contact area of phase 
transformed zone is about three times the uniaxial compressive 
transition stress^^. In addition, with an asymmetrical trend of mar- 
tensite deformation in a uniaxial compression-tension test, the uni- 
axial compressive transition stress is about 1.5 times the tensile 
transition stress^^. 

Eq. 1 is used to define the dissipated energy of NiTi SMA during a 
complete round of the forward and reverse phase transition which is 
defined as^^: 



Hys = H-£rr 



(4) 



where is the maximum phase transition strain. 

The parameters in Eq. 1-4 can be limited to an isothermal loading 
condition (slow loading rates) with constant a am and Gma- In real 
applications, the alloy usually experiences a non-isothermal con- 
dition in the indented bulk during the loading-unloading stages. 
This has a significant effect on the parameters in Eq. 1-4 due to 
the high sensitivity of NiTi SMA to temperature^^. During a cyclic 
loading, the variation of temperature is the result of a change in the 
rate, either loading rate or heating/cooling rate^. As such, the vari- 
ation of hysteresis energy is governed by two main sources: the 
loading-unloading rate in each cycle and the total accumulated heat 
during a cyclic set up. 

In a nanoindentation test with the dominancy of conduction, the 
critical time parameter for heat transfer is defined as^^: 



^critical ~ 



ted 



ti, ted are the characteristic loading time and conduction time defmed as: 

fL=^ (7) 

F 

a' is the state coefficient (depending on the tip and the initial condition), 
X is the heat capacity per unit volume, R is the tip radius, k is the thermal 



conductivity, e is the equivalent indentation strain for a spherical tip 

defined as , where a is the indentation contact radius^^, [i is the rate 

coefficient^^, F^y^^x is the maximum load and F is the loading rate. The 
amount of heat transfer energy also depends on the contact area with 
the tip (heat transfer through the tip) and the interface area (heat 
transfer through the bulk). In a tensile test, the amount of phase trans- 
ition volume and its resultant heat release are limited to a certain 
amount for the tested piece's geometry^^. In contrast, during an indenta- 
tion test, with an unlimited bulk beneath the tip, the time for releasing 
latent heat in Eq. (7) is determined by the maximum load and the 
loading rate (not by g^). 

The critical and characteristic times in Eq. 5-7 determine the 
isothermal, rate dependent and adiabatic performance zones for a 
thermomechanical SMA^^'^^'^^. The critical time (Eq. 5), as one of the 
main parameters in a dynamic loading on thermomechanical mate- 
rials, introduces a balance limit between the conduction time (heat 
transfer time) and the loading time (heat release time). In general, for 
an isothermal condition when tcd^ the temperature remains 
unchanged equal to the initial loading temperature. Preventing the 
temperature increase by conducting most of the released latent heat 
concludes in an unchanged phase transition stress. This results in a 
higher amount of phase transformed volume. A longer cooling time 
(route 1-2-3 in the inset Fig. 2) and a higher thermal relaxation in an 
isothermal condition stabilize the temperature in the transformed 
zone as well as in the non-transformed zone which is subsequently 
transformed during the loading process. The stabilized temperature 
has been reported to be equal to the initial temperature in an iso- 
thermal macroscale tensile condition with a convection domi- 
nancy^\ In the rate-dependent or adiabatic condition, when 
tcd^ tiy the temperature builds up. This phenomenon is due to 
either faster loadings (high amount of heat release and less conducted 
heat) or non-isothermal consecutive cyclic loadings (accumulated 
untransformed/unabsorbed heat). This leads to an increase in the 
temperature and the phase transition stress and, therefore, less phase 
transformed volume and dissipated energy. In the next section, as 
technologies to trace the phase transition and to measure the tem- 
perature of the nanoscale bulk during an indentation test are not yet 
available^^, we performed a numerical simulation at different tem- 
peratures to support the above rationale. 

Numerical simulation-constitutive theory. A phenomenological 
3D constitutive model was implemented into ABAQUS 2010. The 
model is based on the general inelastic frame and the total strain 
tensor £ is disassembled into elastic strain the transformation 
strain £^ and martensite plastic strain £^: 



= £^ + £^'^ + £^ 



(8) 



The elastic strain is assumed to be related to the stress by the effective 
elastic modulus tensor D^^(z): 



(9) 



where <t is the stress and z is the martensite volume fraction. The 
model is constructed from the assumption that the martensitic 
transformation mechanism is associated with an active reorientation 
process, described by the following equation: 



dF{<j,0) 
da 



(10) 



F{g,9) is the loading function that drives the martensitic transfor- 
mation. In this simulation, the Drucker-Prager type of loading 
function is used for the forward and reverse transformation: 



F(a,^) = ||s||+3aP-^AM^ 



(11) 
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Figure 4 | Finite element simulation results for indentations using a 3.61 ^m radius tip: The volume fraction of martensite phase at (a) 20 °C (b) 35 °C 
(c) 50°C, and the plastic strain at (d) 20°C, (e) 35°C, (f) 50°C. 



F((T,^) = ||s||+3aP-^?MA^ (12) 

s = a — PI is the deviatoric stress tensor where P=-g:1 is the 

3 

hydrostatic pressure, I is the second rank unit tensor, and a reflects the 
anisotropic transformation responses in tension and compression 
processes. 

After finishing the forward transformation, the plastic deforma- 
tion occurs in the stress-induced martensite if the applying stress 
reaches the plastic yield stress of the martensite phase. Assuming 
that the plastic deformation obeys the von Mises yielding condition, 
the plastic strain can be shown as: 

d(5 V2^||s|| ^ ^ 

where p is the effective plastic strain, k' is the plastic multipUer, and 
Fy is the martensite plastic yield function. 

The material properties are selected as: same elasticity modulus of 
phases: ^ ^ 31 GPa, same Poisson's ratio of phases: = Vm 
= 0.33, Bam = 6.0 MPa/K, Bma = ^-6 MPa/K, maximum phase 
transition strain: = 0.048, yield stress of the martensite phase: 
0^ = 900 GPa, testing temperatures: 20°C, 35°C, 50°C, a = 0.178, 

g'j^^ = 385 MPa, c/^^ = 402 MPa, = 215 MPa, c/^^ = 198 MPa, 
where superscript s and / show the start and finishing stress of the 
forward transformation (AM) or the reverse transformation (MA). 



In this finite element simulation, the spherical indentation was 
modeled as a 2D axisymmetric problem. Roller boundary condition 
is imposed on the symmetric axis of the sample and the bottom is 
constrained in both radial and axial directions. The minimum ele- 
ment length in the contact zone is 12.5 nm. To avoid the influence of 
boundary condition, the side length of the sample is thirty times 
larger than the radius of the indenter. A preliminary mesh sensitivity 
analysis was performed to ensure that the simulated results were 
insensitive to the element size in the contact zone. The results of 
the above simulation are presented in Fig. 4 (a-f). These figures show 
that the transformation volume and the martensitic plastic zone 
decrease with an increase in the testing temperature. 

Discussion 

The hysteresis energy variation in Fig. 3 can be explained by Eqs. 1-3. 
By raising the temperature both a am W) and a ma W) increase. 
During this stage, experimental reports show that at higher rates 
the growth of Gma W) is higher than the growth of Gam W) due to 
the relationship 0 < Bam < ^ma^^'^°- Thus, comparing two isothermal 
and non-isothermal setups, either at a higher loading-unloading rate 
or as a result of consecutive cycles, H and Hys (Eq. 4) decrease due to 
the higher temperature caused by the accumulative heat and a lower 
cooling time^ In fact, the transfer time for 2000 /nN/s is 10 times 
lower than the isothermal rate 200 fiN/s. This leads to less heat 
transfer after each cycle and, therefore, more total heat accumulation 
after several cycles. While the cycle sequence increases, the phase 
transition stress also increases due to the increase in temperature 
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(Eq. 2). By increasing the phase transition stress and temperature, the 
phase transformed volume decreases consequently (see Figs. 4 (a-c)). 
This phenomenon provides a balance between the phase transition 
volume and the accumulation heat which results in a constant hys- 
teresis energy. In other words, there is a balance between the con- 
duction time and the released time that concludes in an unaltered 
trend of hysteresis energy after several cycles. 

According to Figs. 4 (d-f), when the temperature increases the 
plastic deformation also reduces. This fact can be seen in Fig. 1 where 
during the cyclic loadings the residual depth of each cycle decreases. 
Fig. 4 also shows that the volume of the bulk with 100% martensite 
phase (red color zone) is much larger than the volume of the bulk 
with high plastic deformations, e.g., > 2%. This fact confirms that the 
main determining characters for SMA's behavior during the cyclic 
spherical nanoindentation are phase transition and latent heat and 
not plasticity^^. As such, the dislocation evolution has relatively sim- 
ilar influence in both setups. 

At the isothermal low rate, the hysteresis energy captures the 
stabilized amount after an increasing trend. With the absence of 
temperature increase, the hysteresis increase is explained by the 
reorientation phenomenon in favour of the phase transition in sub- 
sequent cycles^°'^^. For the isothermal loading rate, the phase trans- 
ition stress is lower (more achievable by a similar applied load, 
10000 fiN) than the transition stress level for the non-isothermal 
loading rate. Therefore, after stabilization the final hysteresis energy 
for the rate 200 fiN/s is expected to be higher than the hysteresis 
energy for the rate 2000 jdN/s. In brief, with a similar initial con- 
dition, the atomic and microstructure degradations influence the 
cyclic-loading results for the isothermal condition (Fig. 2) and the 
two time scales govern the hysteresis trend for the higher rate setup 
(Fig. 1). 

In comparison to the lower rate data. Fig. 3 shows a higher amount 
of hysteresis energy in the first few cycles for the higher rate setup. 
There are three possible reasons for this. First, with an increase in the 
loading rate, stress-induced phase transformed nucleation occurs in 
an increasing number of different locations (Figs. 5 a, b). 

Secondly, the phase transformed volume is propagated by an 
interface with the velocity of v^^: 



rph 



(14) 



Vph represents the maximum radius of the bulk underneath the tip 
where the stress is higher than the nucleation phase transition stress. 
Eq. (14) implies that for a similar size of the loaded bulk, faster load 
(lower tl) has a higher velocity of interface propagation. This con- 
cludes in a higher phase transformed volume with the resultant 
larger hysteresis area at higher loading rates for the first few cycles 
in which the heat has not been significantly accumulated by the 
following cycles. While the heat is accumulated, the velocity 
decreases subsequently. 

Fig. 5 schematically provides a hypothesis for the nucleation and 
propagation of the phase transition volume in a nanoindentation 
process. When the loading rate is in an isothermal condition, the 



nucleation and propagation of the phase transition volume mostly 
occur with one or two interfaces at a low speed (Fig. 5a)^\ As the 
loading rate increases, the number of nucleation sites and the inter- 
face velocity increase (Fig. 5b). The domain spacing between the 
phase transition sites and the number of nucleation sites have been 
shown to be proportional to and £ in the dominancy of the 
convection (tensile testing) Similar nucleation and propagation 
sites in Fig. 5b were recently shown as local micro-plastic hardening 
and twinning by a MD simulation for indentation experiments^^. In 
the present study, compared to the smooth indentation curves for the 
isothermal rate setup (Fig. 2), the jagged shape of the indentation 
curves (Fig. 1) for the non-isothermal rate setup is due to numerous 
new nucleation sites beneath the tip while the load is applying. These 
sites relax the stress and increase the strain at the time of nucleation. 

It should be noted that the number of nucleation sites also depends 
on the defect zones such as inclusions, existent dislocations and 
precipitations where the stress concentration nucleates the phase 
transformation. The latter two are shown in Fig. 6 using AsB and 
TEM. This figure indicates two major single crystal precipitations 
located at different places of the texture. The dislocation congestion 
around the precipitations with a hard phase is clearly seen in the 
TEM image in inset Fig. 6. The existence of these structures in the 
stressed bulk can be one of the rationale of having either irregular or 
dissimilar jagged events (nucleated phase transition volume) 
throughout an indentation curve in Fig. 1. 

Furthermore, the activation volume in a thermally activated 
mechanism has the following relationship^^: 



V* = V3k0 



d\n8\ y/2>ke 



(15) 



where £, m and k are strain rate, rate sensitivity and Boltzman's 
constant. In this equation, the phase transition stress is assigned as 
(7 for NiTi SMA for the applied stress using a spherical tip. By Eq. 15 
the previous rationale, higher hysteresis energy in the first few cycles, 
is validated as the activation volume is larger at higher rates. This can 
be associated with the higher initial propagation velocity of interface 
and the larger radius of stressed bulk (Eq. 14). 

Thirdly, using Eq. 2 and 3, on increasing the temperature, Eq. 15 
degenerates into the following equation: 




dislocation congestion 



□ 




g with the velocity v 



Figure 5 | Nucleation sites of the phase transition volume and their 
propagation under the cyclic indentation loading with a lower rate (a), 
and with a higher rate (b). 



Figure 6 | TEM and AsB images of precipitations and the presence of 
dislocations around them. 



SCIENTIFIC report; | 3 : 341 2 | DOI: 1 0.1 038/srep0341 2 



5 



AV* = 



V3k 
mb 



(16) 



Ay* represents either the activation volume variation after each 
consecutive cycle or the total change in the activation volume after 
10 cycles. Eq. 16 shows that unlike ordinary materials with no phase 
transition, if the temperature increases during the cyclic loading, the 

activation volume decreases with the ratio - (<1) where b is the 

b 

Clausius-Claperon coefficient shown in Eqs. 2 and 3. This is an 
additional rationale for a decreasing hysteresis during the cyclic 
loading with a higher rate (Fig. 3). During the cyclic loading with 
an isothermal rate, i.e., 200 jiN/s, the activation volume remains 
unchanged due to the constant temperature. 

It is not coincidental that accumulated heat during the cyclic load- 
ing has a close correlation with the initial condition, e.g., initial 
temperature and hysteresis energy. In cyclic tensile loadings, the 
stable mean temperature and its fluctuation depend on the initial 
temperature, its initial hysteresis energy and the characteristic time 
scale of convection and the loading-unloading time. The total trans- 
ferred heat to the ambient remain unchanged after several cycles. 
Similarly, in nanoindentation, the final condition of each cycle, 
which is the initial condition for the subsequent cycle, plays an 
important role in raising the accumulated heat at the end of 10 cycles. 

In summary, the behavior of NiTi SMA was studied under two 
different cyclic loading setups and a hypothesis was provided on the 
nucleation of the phase transformation volume in a nanoscale 
stressed bulk. To obtain an improved cooling process, a cycUc load- 
ing with a lower rate and a holding time was applied on NiTi SMA. In 
this setup, the hysteresis energy increases from the initial cycles to a 
stabilized magnitude in subsequent cycles. This trend is explained by 
the reorientation that situates the alloy in a favourable modal energy 
for the subsequent cycles. In the second cyclic loading setup with a 
higher rate and a weaker cooling process, the time balance between 
the heat release and the heat conduction determines the trend of 
lowering the hysteresis energy in subsequent cycles. 

Methods 

Using a differential scanning calorimeter (DSC, TA-QIOOO), the austenite finish 
temperature (Af) of NiTi SMA was measured as 19°C. The material was initially in the 
austenite phase and the martensite phase was induced by the applied compressive 
stress with a superelastic behavior. The nominal alloy composition as determined by a 
high resolution X-ray diffraction spectroscopy (XRD, PW1825) was 56.4% Ni and 
43.6% Ti by weight. The grain size was in the range of 50-100 nm as observed by a 
transmission electron microscopy (TEM, 200 kV: JEM-2100F). Angle selective 
Backscattered electron (AsB) images of orientation and precipitation were obtained 
using a Zeiss Supra 55VP. The sheet was cut into pieces (5 mm X 5 mm) and polished 
using a series of silicon carbide and diamond papers until the average surface 
roughness was less than 6 nm measured using a 3D surface profiler (SPM, NT3300). 
In order to achieve a well- controlled compression process, an in situ nanoindentation 
machine, a Hysitron nanoindentor (TI 900 Triboindenter) was used. Before carrying 
out tests, the instrument was carefully calibrated using standard procedures. The 
nanoindentation tests were conducted at room temperature (23°C) in a quasi-static 
mode using a spherical diamond tip with a 3.61 |im radius. A constant maximum 
load of 10000 |iN was applied in all tests. The indentation load was controlled by an 
ultrafine piezoelectric beam in the machine. 
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